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We analyze a theoretical model for energy and electron transfer in an artificial photosynthetic 
system. The photosystem consists of a molecular triad (i.e., with a donor, a photosensitive unit, and 
an acceptor) coupled to four accessory light- harvesting antennas pigments. The excitation energy 
transfer from the antennas to the artificial reaction center (the molecular triad) is here described 
by the Forster mechanism. We consider two different kinds of arrangements of the accessory light- 
harvesting pigments around the reaction center. The first arrangement allows direct excitation 
transfer to the reaction center from all the surrounding pigments. The second configuration transmits 
energy via a cascade mechanism along a chain of light-harvesting chromophores, where only one 
chromophore is connected to the reaction center. At first sight, it would appear that the star- 
shaped configuration, with all the antennas directly coupled to the photosensitive center, would be 
more efficient. However, we show that the artificial photosynthetic system using the cascade energy 
transfer absorbs photons in a broader wavelength range and converts their energy into electricity with 
a higher efficiency than the system based on direct couplings between all the antenna chromophores 
and the reaction center. 

PACS numbers: 



I. INTRODUCTION 

The reaction centers of natural photosystems are sur- 
rounded by a number of accessory light-harvesting com- 
plexes pp. These light-harvesting antennas absorb sun- 
light photons and deliver their excitation energy to the 
reaction center, which creates a charge-separated state. 
The photosystem of green plants is made up of six pho- 
tosynthetic accessory pigments: carotene, xanthophyll, 
phaeophytin a, phaeophytin b, chlorophyll a, and chloro- 
phyll b pp. Each pigment absorbs light in a differ- 
ent range of the solar spectrum. As a result, the an- 
tenna complex significantly increases the effective fre- 
quency range for the light absorption, resulting in a 
highly-efficient photocurrent generation. In the presence 
of excessive sunlight the antenna complex can reversibly 
switch to the photo-protected mode, where harmful light 
energy is dissipated. 

The efficient performance of natural photosystems mo- 
tivates researchers to mimic their functions by creating 
photosynthetic units, which are combined to antenna 
complexes with artificial reaction centers. For exam- 
ple, a light-harvesting array of metalated porphyrins has 
been developed in Ref. [2J. This array absorbs light and 
rapidly transfers the excitation energy to the reaction 
center, so that the porphyrin (P)-fullerene (Ceo) charge- 
separated state, P + -C^g, is formed with a quantum yield 
~ 70%. Mixed self-assembled monolayers of the ferrocene 
(Fc)-porphyrin-fullerene molecular triad and the boron 
dipyrrin (B) dye have been made in Ref. 0] with the 
goal to examine both energy and electron transfers in 
the artificial reaction center (Fc-P-Cgo), coupled to the 
light-harvesting molecule B. A quantum yield of ~ 50% 
for photocurrent production at the wavelength 510 nm 
and a quantum yield of ~ 21% at the wavelength 430 nm 
have been reported [3]. 



Recently, a more efficient, sophisticated and rigid 
antenna-reaction system has been designed in Refs. [3[6]. 
This system includes three kinds of light- absorbing chro- 
mophores: (i) bis (phenylethynyl) anthracene (BPEA), 
which absorbs around 450 nm wavelength (blue region); 
(ii) borondipyrromethene (BDPY), having a strong ab- 
sorption at 513 nm (green region); and (iii) zinc tetraaryl- 
porphyrin, which absorbs both at 418 nm and at 598 nm. 
This study reports ~100% quantum yield for the excita- 
tion transfer and ^95 % quantum yield for the generation 
of the charge-separated state P + -Cg . 

Energy transfer mechanisms in the multi-chromophoric 
light-harvesting complexes of bacterial photosystems 
(e.g., excitation-transfer between chlorophyll molecules 
in the Fenna-Matthews-Olson (FMO) protein complex) 
have been studied elsewhere (see, e.g., Ref. [7HlQ]). 
Those works [TlU0| have mainly focused on the quantum 
effects in excitation-transfer across the bacteriochloro- 
phyll units. 

Theoretical studies of antenna-reaction center com- 
plexes can be useful for a better understanding of, and 
for optimizing light-to-electricity conversion, as well as 
for developing new and efficient designs of solar cells. 
Recently we have theoretically analyzed pj] the light-to- 
electricity energy conversion in a molecular triad (Fc-P- 
C6o) electronically coupled to conducting leads. It was 
shown that the Fc-P-Ceo triad can transform light en- 
ergy into electricity with a power-conversion efficiency of 
order of 40%, provided that the connection of the triad to 
the leads is strong enough. It should be noted, however, 
that this prototype solar cell absorbs photons with ener- 
gies in close proximity to the resonant transition of the 
central porphyrin molecule. Therefore, a major fraction 
of the sunlight spectrum is not converted to the electrical 
form by this device. 

In this paper we examine a theoretical model for 
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FIG. 1: (Color online) Schematic diagram of an artificial pho- 
tosystem comprised of a molecular triad (D-P-A) and four 
additional light-harvesting complexes (Anl, An2, An3, An4). 
Here, D = donor, A = acceptor, and P = photosensitive part 
(porphyrin). The molecular triad D-P-A is inserted between 
two electrodes (leads) L and R. Energy exchange processes 
are denoted by straight red arrows. The green bent curved 
arrows describe electron pathways, L — > D — s> P — >A— > 
R, via the molecular triad, (a) The photosensitive part, P, 
of the molecular triad is surrounded by four accessory light- 
harvesting complexes Anl, An2, An3, and An4. In this case 
the surrounding antenna complexes can transfer excitations 
to the reaction center directly, (b) Here the antenna com- 
plexes form a linear chain coupled to the reaction center via 
nearest-neighbor couplings. 



the light-to-electricity energy conversion by a molecu- 
lar triad, which is surrounded by four additional light- 
harvesting antenna complexes. We show that this ar- 
tificial photosystem is able to generate a photocurrent 
with a quantum yield of the order of 90 % (when the re- 
organization energy for the Forster transfer is relatively 
high) absorbing photons in a wide range (420-670 nm) 
of the solar spectrum. We consider two different con- 
figurations for the antenna complexes: (a) where each 
light-harvesting molecule is independently connected (by 
the Forster energy-transfer mechanism) to the central 
porphyrin (P) molecule of the triad (Fig. la); and (b) 
where the light-harvesting molecules are arranged in a 
line (Fig. lb) , with only one molecule directly coupled to 
the porphyrin and with other molecules forming a chain 
where the energy propagates in a cascade-manner. 

Let us consider a (D-P-A) reaction center surrounded 
by several (say, four) accessory light-harvesting anten- 
nas. Which is the best way to arrange in space these 
accessory antennas? In other words, which network or 
topology would provide more energy from the surround- 
ing accessory antennas to the central reaction center? A 
star-shaped configuration? (with each accessory antenna 
directly coupled to the central reaction center). Or in a 
somewhat opposite configuration: as a linear chain with 
nearest-neighboring couplings between the antennas and 
only one of these coupled to the reaction center? These 
two cases are shown in Fig. 1. 

In principle, very many possible topologies could be 
considered. However, to simplify this analysis, we will 
now focus on two extreme cases, with somewhat opposite 
topologies or networks: a well-connected reaction center 
(directly connected to all four accessory light harvesting 
antennas), and the opposite case where the central re- 
action center is coupled to only one antenna, which is 
now part of a linear chain. These two extreme-opposite 
topologies or networks can be denoted as star-shaped 
(shown in Fig. la) and linear-chain (Fig. lb), respec- 
tively. 

At first sight, it would seem that the star-shaped con- 
figuration, with each antenna directly connected to the 
reaction center, can provide far more energy to the reac- 
tion center, due to the multiple connections between the 
central unit, and the surrounding antennas. However, 
an energy bottleneck can arise even in this apparently- 
optimal topological configuration. Each antenna is here 
assumed to operate optimally in a limited, and shifted, 
frequency range. In other words, we are not assuming all 
accessory light-harvesting antennas to be equal to each 
other. In general, various antennas can operate in dif- 
ferent frequency ranges, and this difference is crucial in 
our analysis. Recall that the photosystem of green plants 
is made of six light-harvesting accessory pigments: each 
pigment absorbing light in a different range of the solar 
spectrum. Thus, in the star-shaped topology shown in 
Fig. la, antennas with a large energy mismatch would 
not provide energy to the central reaction center. Only 
the surrounding accessory antennas that have an approx- 



3 



imate energy match to the central reaction center would 
transfer energy. This energy bottleneck works against 
the star-shaped network shown in (a). 

Thus, the two issues considered here are: (1) how to 
physically arrange antennas around the central reaction 
center, and also (2) how to arrange these in energy-space, 
so to speak. The first issue is topological and focuses on 
the network connectivity in real space: for instance, how 
many accessory antennas are connected to a central reac- 
tion center. The second issue refers to the energy match 
(or mismatch) between neighboring antennas, and be- 
tween them and the central reaction center. A large en- 
ergy mismatch between any connected units in the chain 
would preclude energy transfer between them. This ap- 
proximate "energy matching" issue between neighboring 
units is equally important to keep in mind, not just the 
real-space topological arrangement of the units. 

As mentioned above, Fig. 1 shows the connectivity be- 
tween the different units: star-shaped topology in (a), 
and linear-chain configuration in (b). Moreover, the col- 
ors there represent, very schematically, the energy range 
where each unit operates optimally. The linear chain 
shown there has antennas arranged in a way that nearby 
units operate in approximately similar energy ranges. 
This energy-matching issue perhaps is not very clear in 
Fig. 1, even when seen in color. The energy scales are 
shown far more explicitly in Fig. 2. Figure 2b clearly 
shows that the linear chain model considered here op- 
erates via an energy cascade, or linear chain-reaction. 
Like a line of falling dominoes, one event triggering the 
next one, in a sequential manner, with small energy mis- 
matches between successive energy transfer events. As 
shown in Fig. 2b, the more energetic antenna is located 
far away from the reaction center, and it is coupled to an 
antenna with a slightly lower energy, which is coupled to 
another antenna with an even slightly lower energy, and 
so on, moving energetically "downhill" along the chain. 
Thus, the neighboring antennas must be so both in real 
space, and also in "energy space", to allow for the effi- 
cient transfer of energy between them. Thus, "proxim- 
ity" between units must be in two spaces: real space and 
energy space. 

This article is organized as follows: in Section II we 
outline a model for the artificial reaction center (molec- 
ular triad) coupled to the antenna complex. We briefly 
describe our mathematical methods in Sec. III. The pa- 
rameters are listed in Sec. IV. In Sec. V, we numerically 
solve the master equations and analyze the energy trans- 
fer process. Conclusions are presented in Sec. VI. The 
methods used are described in more details in the online 
supplementary material. 



II. MODEL 
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FIG. 2: (Color online) Energy diagram of the antenna com- 
plexes Anl, . . ., An4, energetically coupled to the reaction 
center (RC), D-P-A, for the case (a) where each light- 
harvesting molecule is directly connected to the porphyrin 
molecule, P, of the RC; and (b) where the energy transfer 
occurs via a linear-chain of antenna with nearest-neighbor 
energy exchange. Note that the four antenna chromophores 
in (b) are arranged in a way to have a relatively small energy 
difference between neighboring units. 



We start with a schematic description of the energy 
and electron transfer processes in an artificial reaction 
center D-P-A (Donor-Porphyrin-Acceptor) combined 
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with four antenna chromophores: Anl, An2, An3, An4 
(see Figs, la and lb, showing two different configurations 
for these antennas: star-shaped configuration in la and 
chain in lb). The photosensitive molecular triad, D-P- 
A, is inserted between two electron reservoirs (electrodes) 
L and R. The donor, D, is coupled to the left lead, L, and 
the acceptor, A, is connected to the right lead, R. As in 
Refs. [21 [TT], the donor and acceptor molecules (e.g., fer- 
rocene and fullerene), are connected to each other via the 
photosensitive molecule (porphyrin, P). This molecule is 
surrounded by four Accessory Light-Harvesting Pigments 
(ALHP). Figure la corresponds to the situation where 
all pigments are directly coupled to the photosensitive 
part (P) of the molecular triad. In this "star-shaped" 
geometric arrangement, all the ALHP can directly trans- 
fer excitations to the photosensitive part of the molecu- 
lar triad. Figure lb corresponds to the case where the 
light-harvesting pigments form a chain, which transfers 
energy where the excitation moves from one pigment to 
the next one via nearest-neighbor couplings: An4 — > An3 
— > An2 — > Anl — > P. This cascade-like excitation trans- 
fer occurs in an energetically-downhill direction, akin a 
one-dimensional chain reaction or domino- effect. 

Figures 2a and 2b present the energy diagrams of the 
photosystems described in Figs, la and lb, respectively. 
The electron transfer chain, L^D— >P^ P* — >A— >R, is 
the same for both configurations (a) and (b), and both 
begin on the left lead, L. The electrochemical potentials 
of the left (L) and right (R) electron reservoirs are de- 
termined by the parameters fip and fip, with > /il- 
Since the energy level Ep, of the donor D is lower than 
the potential /ip of the left lead, Ep < fip, electrons 
can move from the L-reservoir to the level D and, after- 
wards, to the low- lying ground energy level, Ep, of the 
porphyrin. When absorbing a photon, the electron in 
the porphyrin molecule jumps from its ground state P 
to its excited state P*. A subsequent electron transfer 
from P* to the acceptor A is driven by a negative en- 
ergy gradient, (Ea — Ep*) < 0. In view of the relation: 
Ea > Hr, the electron in A is finally transferred to the 
right, R, electron reservoir. This is the light-induced elec- 
tron transition in the porphyrin molecule, which results 
in an energetically-uphill electron flow in both photosyn- 
thetic systems (a) and (b). 

Even though these systems (a) and (b) have the simi- 
lar electron-transport chains, their light-harvesting com- 
plexes are arranged quite differently. Each of these com- 
plexes, An = Anl, . . ., An4, can be characterized by a 
ground, Ea u , and an excited, Ea u *, energy levels with 
an energy difference uiAn = Eau' — Eau- Hereafter, 
we assume that Ti = 1 and kp = 1- For the light- 
harvesting complex (a) (see Figs, la, 2a) all frequencies 
oJAnii ■ ■ • i^Ani should exceed the porphyrin transition 
frequency, tup = Ep*—Ep. In this case the energy of pho- 
tons collected by each individual antenna can be trans- 
ferred directly to the photosensitive part of the artificial 
reaction center (porphyrin molecule). However, in the 
light-harvesting complex (b) (see Figs, lb and 2b) only 



the antenna Anl is coupled (by a Forster mechanism) 
to the porphyrin, whereas the other light-harvesting pig- 
ments form a linear chain that transfers energy down- 
hill along the chain: An4 — > An3 — > An2 — > Anl — > 
P. This energy transfer can be energetically-allowed pro- 
vided that LOAni > ^An3 > W An2 > ^ Aril > Up. In SeC. V 

we compare these two artificial photosystems and deter- 
mine which arrangement of the antenna complexes pro- 
vides more energy to the reaction center. 

III. METHODS 

The electron flow through a molecular triad coupled 
to two electron reservoirs can be described with meth- 
ods of quantum transport theory and the theory of open 
quantum systems [T^HTS] . In addition to four sites (D, 
P, P*, A), describing the molecular triad, we introduce 
four pairs (Anl, Anl*, . . ., An4, An4*), which character- 
ize the ground and excited states of the light-harvesting 
antennas. The whole system, the molecular triad plus 
four antenna complexes, can be analyzed within a math- 
ematical formalism presented in Supporting Information 
and also in Ref. [TT] . 

The total Hamiltonian of the system includes the fol- 
lowing components: 

(i) energies of the electron sites and leads, as well as the 
Coulomb interactions between the electrons located on 
different sites of the triad; 

(ii) tunneling couplings between the electron sites on the 
triad and the electron reservoirs; 

(iii) electron tunneling between the electron sites belong- 
ing to the molecular triad; 

(iv) coupling of electron sites and antenna complexes to 
an environment; 

(v) interactions of the porphyrin molecule and antenna 
complexes with an external electromagnetic field (laser 
field) and with 

(vi) blackbody radiation and Ohmic bath, responsible for 
the quenching (energy loss) of the porphyrin and antenna 
excited states. 

(vii) In the case of the design in Fig. la the Hamiltonian 
includes the direct Forster coupling between the por- 
phyrin molecule, P, and the light-harvesting complexes 
Anl, . . ., An4. For the linear-chain configuration shown 
in Fig. lb, the Forster mechanism provides the energy 
transfer between the nearest-neighbors in the antenna 
chain, as well as between the complex Anl and the por- 
phyrin molecule. 

IV. PARAMETERS 

Energy levels and electrochemical potentials: The en- 
ergy levels of the Fc-P-Ceo molecular triad are Ep, = 
-510 meV, E P = -1150 meV, E P * = 750 meV and 
Ea = 620 meV. These values are obtained by esti- 
mating the reduction potentials (using a reference elec- 



5 



trode Ag/AgCl) of ferrocene (D), porphyrin (P, P*) and 
fullerene (A) molecules [16]. For the electrochemical 
potentials of the left (/ip) and the right ([Ir) leads, 
we choose the following values: [ip = —410 meV and 
Hp — 520 meV, with the electrochemical gradient A/x = 
Hr. — Hl = 930 meV. 

Coulomb interactions: The spatial separations be- 
tween D-P, P-A and D-A are of order of 1.62 nm, 1.8 
nm and 3.42 nm, respectively [16 . The Coulomb energies 
mop, uda and up a can be calculated with the formula 
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where, {ij} = {DP}, {DA}, {PA}, and e is the vac- 
uum dielectric constant. For e ~ 4.4, the Coulomb inter- 
action energies are up, p = 200 meV, u da = 95 meV and 
up a — 180 meV. 

The Forster coupling, Vp, between the photosensitive 
molecules I and V is proportional to the product of the 
dipole moments of these molecules, er\ and erv, and in- 
versely proportional to the cubic power of the distance 
(Rw ) between them [17] : 
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For the case where ~ r; ~ 0.3 nm, R ~1 nm, and at 
e ~ 4, the Forster coupling is about Vp ~ 65 meV. 

Tunneling amplitudes: We have assumed the 
ferrocene-porphyrin and porphyrin-fullerene tunneling 
amplitudes are about ~ 3 meV, so that 
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For the tunneling rates between the left lead and fer- 
rocene (r^) and between the right lead and fullerene (r^) 
we choose the following values [11]: Tp/h = 1800 /is -1 
and T R /h = 180 ^s" 1 . 

Radiation leakage and quenching rates: We take 
the following estimates for the radiation leakage time: 
t p*->p = T P'^>-D — T A-fP ~ 0.4 ns. Similar estimates 
have been used for the radiation leakage timescales of 
the antenna molecules. For the quenching (or energy- 
loss) time of the porphyrin excited state P* we use the 
value: r qucn ~ 0.1 ns. 

Reorganization energies for the electron and energy 
transfers: For the molecular triad analyzed in Ref. [TT] 
we obtain the relatively high power-conversion efficiency, 
•q ~ 42%, provided that the donor-porphyrin and 
acceptor-porphyrin electron transfer reorganization en- 
ergies are about Apjp ~ 600 meV and Aap ~ 100 — 
400 meV. A much smaller value, App» ~ 100 meV, is as- 
sumed for the light-induced electron transitions between 
the ground (P) and excited (P*) levels of the porphyrin 
molecule. The Forster energy transfer between the light- 
harvesting molecules and between these molecules and 
the porphyrin is also accompanied by an environment- 
reorganization process, which can be characterized by a 



smaller energy scale, Ap < 100 meV (see, e.g., the energy 
transfer in the B850 complex [18], where Ap is assumed 
to be about 30 meV). 

Hereafter, we assume that the external light source 
has a fixed intensity, I = 100 mW/cm 2 , and that the 
environment is kept at the room temperature, T — 298 K. 
We also assume that the reorganization energy for the 
Forster energy transfer, Ap, is about 100 meV, unless 
otherwise specified. 

We analyze the Fc-P-Ceo molecular triad, where the 
ferrocene molecule, Fc, is attached to the gold surface 
(left lead, L), and the fullerene, Cqo, is in contact with 
an electrolyte solution (right lead, R) filled with oxygen 
molecules, which are able to accept electrons from the 
Ceo molecules. 



V. RESULTS AND DISCUSSION 

We derive and solve numerically a set of master equa- 
tions for the probabilities to find the system in a definite 
eigenstate of the basis Hamiltonian. This is explained in 
the online supplementary material. After that, we calcu- 
late: (i) the energetically-uphill electron current through 
the triad, (ii) the energy of the photons absorbed by the 
triad and by the light-harvesting molecules. This allows 
us to determine a quantum yield and power-conversion 
efficiency of the system (see all definitions in the support- 
ing information). 



A. Photocurrent through the molecular triad 
directly coupled to four porphyrin light-harvesting 
molecules 



Here we consider the situation where both the reaction 
center and the antenna complexes are made of porphyrin 
molecules with the geometrical arrangement shown in 
Fig. la. This arrangement allows direct energy transfer 
from each light-harvesting chromophore to the reaction 
center. 

In Fig. 3 we plot the photocurrent through the triad as 
a function of the wavelength of light for different values 
of the Forster coupling strength : Vp = 0, 0.1, 1, 10 (in 
meV) and for the above-mentioned set of parameters of 
the system. It is apparent from Fig. 3 that the magnitude 
of the light-induced pumping current at A = 620 nm is 
significantly enhanced (about 5 times larger when Vp = 
10 meV) by the antenna system. However, the spectral 
range of the light absorption remains the same as for the 
detached porphyrin reaction center (see Fig. 3, where the 
red curve with square symbols describes the photocur- 
rent through the molecular triad completely disconnected 
from the antenna chromophores, Vp = 0). We also find 
that the quantum yield, $, taken in the middle of the res- 
onant peak (A = 620 nm), non-monotonically depends 
on the Forster coupling strength Vp measured here in 
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FIG. 3: (Color online) Photoinduced electron current I e ( 
number of electrons, in one ms, pumped from the L to the 
R lead) versus the wavelength of the incident light for a pho- 
tosystem with four antenna complexes, which are made of 
porphyrin molecules. Both the antenna complexes and the 
reaction center absorb at the same wavelength. The whole 
complex now absorbs more photons than the single porphyrin 
molecule; thus, pumping many more electrons from the left 
(with fiL = —410 meV) to the right (with [ir = 520 meV) 
electron reservoir. The other parameters are listed in the text 
(see Sec. IV). The peak in the current increases for larger val- 
ues of the Forster coupling strength Vf- We also studied (not 
shown here) higher values of Vf, but these produced the same 
results as Vf = 10 meV. Thus, this value of Vf provides a 
saturation in the electron current. The resonant peak here is 
A = 620 nm. 
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FIG. 4: (Color online) Electron current I e and quantum yield 
$ as functions of the wavelength A of the external radia- 
tion for the configuration shown in Fig. la, where two BPEA 
and two BDPY antenna chromophores are directly coupled 
to the centrally-located reaction center. The Forster cou- 
pling constant, Vf, which is assumed to be the same for ev- 
ery chromophore-RC connection, takes four values (in meV): 
Vf — 0, 1, 25, 50. For other parameters see Sec. IV. Note 
that the electron current and the quantum yield grow for in- 
creasing values of the Forster coupling energy strength Vf- 
More importantly, the direct coupling (Fig. la) suppresses 
the peak at A= 450 nm (Fig. 4a), which is present in the 
linear chain configuration (Fig. lb), as shown in Fig. 5a 



Now we consider a different case: an antenna system 
comprised of two BDPY and two BPEA molecules. The 
BDPY molecule has the maximum absorbance in the 
green region (at 513 nm) of the solar spectrum, where 
neither BPEA (with maximum absorbance at 450 nm) 
nor the porphyrin, which absorbs at 620 nm, have max- 
ima of absorption spectra. It should be noted that a mul- 
tichromophoric hexad antenna system having three light- 
absorbing BDPY, BPEA and porphyrin has been devel- 
oped in Ref. [5] . The generation of the charge-separated 
state, P + -Cg , with almost 95% quantum yield [5]. In 
our case, the porphyrin unit of the molecular triad is cou- 
pled to the four antenna chromophores (two BDPY and 
two BPEA). 

We consider two situations: (a) where the antenna 
chromophores are directly coupled to the porphyrin unit 



of the reaction center (Fig. la); and (b) where the an- 
tenna chromophores are arranged in line: BDPY — > 
BDPY -s- BPEA BPEA RC, with the nearest- 
neighbor coupling between chromophores. See Fig. lb. 
Thus, the configuration (a) might appear to be energet- 
ically more efficient than (b). However, our calculations 
below indicate that this is not the case, 
(a) For the case of direct connection between the four 
antennas chromophores and the triad (see Fig. la and 
Fig. 4) we calculate a photocurrent and a quantum 
yield, $, as functions of the wavelength of light, at 
Ap = 100 meV, and at five values of the Forster cou- 
pling: Vf — 0, 1, 25, 50 meV. The wavelength depen- 
dence of the current has two maxima centered at 513 
nm and 620 nm. The BPEA molecules, which absorb 
at 450 nm, give a negligible contribution to the cur- 
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FIG. 5: (Color online) Electron current I e and a quantum 
yield $ versus the wavelength of light, A, for the configuration 
shown in Fig. lb, where the excitation energy moves along 
the following chain of light-harvesting molecules: BPEA — > 
BPEA -¥ BDPY -¥ BDPY -> RC. The parameters used here 
are the same as in Fig. 4. However, the current peak at A ~ 
450 nm is present in (a) here, but absent in Fig. 4(a), which 
used a direct-coupling configuration to the central reaction 
center. 



rent since their spectral maxima are too far from the 
absorbance maximum of the porphyrin spectrum. As a 
consequence, the BPEA-porphyrin energy transfer is sig- 
nificantly suppressed at moderate values of the Forster 
reorganization energy, Ap < 100 meV in the range of the 
coupling constants Vf < 50 meV. It follows from Fig. 4 
(see a peak at A = 513 nm) that the BDPY molecules 
start working as efficient light-harvesters only at suffi- 
ciently strong Forster coupling, Vf > 10 meV, to the 
porphyrin unit of the molecular triad. We also note that 
when A ~ 513 nm, both the photoinduced current and 
the quantum yield grow with increasing the Forster cou- 
pling strength, so that the quantum yield, $, can be 
around 48%. In the range of porphyrin absorption (at 
A ~ 620 nm and Vf = 0) the quantum yield is of the 
order of 90%. 

(b) A much broader light spectrum can be converted into 
electrical current in the linear configuration in Fig. lb, 
where the light-harvesting chromophores are arranged 
along a line: BPEA -> BPEA -> BDPY -> BDPY -s- 
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FIG. 6: (Color online) Quantum yield as a function of reor- 
ganization energy, Ap, for the linear-chain nearest-neighbor 
coupling between chromophores, BPEA — > BPEA — > BDPY 
->■ BDPY -> RC, for a the wavelength of light A = 450 nm 
(blue peak in Fig. 5a). This figure is plotted for three different 
temperatures: (a) Very low T = 77 K, (b) room temperature 
T — 298 K and (c) very high T = 500 K as well as for three 
values of the Forster constant Vf = 1, 25, and 50 meV. As 
shown in the figures (b) and (c), increasing the reorganization 
energy Af, can sharply increase the quantum yield. The ex- 
treme low-temperature case in (a) is just a limit case, shown 
for comparison with the higher-temperature cases in (b) and 



RC, with the only one BDPY molecule directly coupled 
to the porphyrin unit of our artificial reaction center (RC) 
(see Fig. lb and Fig. 5). This system is able to collect 
photons in the range of wavelength from 420 nm up to 
650 nm covering a significant part of the visible sunlight 
spectrum. The chain of BPEA and BDPY molecules cre- 
ates an efficient channel, which gradually transmits en- 
ergy from the collectors of high-energy photons (BPEA 
molecules), via the intermediate BDPY antennas, to the 
molecular triad. In Figs. 5a, 5b we plot the photoinduced 
current and the quantum yield versus the wavelength of 
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FIG. 7: (Color online) Quantum yield versus reorganization 
energy, Af, for the green peak (A = 513 nm) of the spec- 
trum (see Fig. 5a) for three values of the Forster constant, 
Vf = 1, 25, 50 meV, and at three different temperatures: 
(a) T = 77 K, (b) T — 298 K and (c) T = 500 K. The other 
parameters are listed in Sec. IV. Figures 5, 6, and 7 focus on 
the linear chain configuration (Fig. lb) with nearest-neighbor 
couplings between chromophores. Figures 6 and 7 show the 
same quantities, but centered at different peaks (A ~ 450 nm 
versus A ~ 513 nm). 



the external radiation, at = 100 meV, and at four val- 
ues of the Forster constants (in meV) Vp = 0, 1, 25, 50. 
For large value of the Forster coupling strength Vp the 
quantum yield <E> reaches ~ 48% in Fig. 4(a) around A = 
513 nm [for the star-shaped topology in Fig. la] and ~ 
30% in Fig. 5(b) [for the linear chain case in Fig lb]. The 
quantum yield in Fig. 5b is lower than the one in Fig. 4a, 
but extends over a wider range of wavelengths including 
the peaks around A ~ 450 nm and A ~ 513 nm. 

It follows from Fig. 4 and Fig. 5, that the configu- 
ration using the chain-like nearest-neighbor coupling be- 
tween light-harvesting chromophores [Fig. lb, design (b)] 
converts much more blue (A = 450 nm) light into elec- 
tricity with a higher quantum yield than the star-shaped 
configuration with direct coupling between the antenna 



chromophores and the RC [Fig. la, design (a)]. In this 
star-shaped configuration, the BPEA molecules, which 
collects blue photons, are not able to transfer their en- 
ergy to the photosensitive unit of the triad due to a sig- 
nificant difference between the energies of the BPEA an- 
tennas (A = 450 nm) and the porphyrin-based reaction 
center (A = 620 nm). 

The antenna-RC energy transfer is facilitated by the 
strong coupling to the environment (when the energy 
difference is high), characterized by the reorganization 
energy Ap, as well as by a tight Forster binding between 
chromophores, which is described by the constant Vf- In 
Figs. 6 and 7 we plot the quantum yield, as a function 
of the reorganization energy, A^?, for three different tem- 
peratures (in K): T — 77, 298, 500, and for three values of 
the coupling constant Vf = 1, 25, and 50 meV. Figure 6 
is related to the blue peak of the absorption spectrum 
(A = 450 nm), whereas figures 7 describes the behavior 
of the green peak (A = 513 nm). The peak centered at 
A = 620 nm is produced by the porphyrin molecule, be- 
longing to the triad, and, therefore, shows no dependence 
on Ap and Vf- 

The Marcus rates, describing the energy transmission 
between the photosensitive elements of the system, de- 
pend (i) on the energy difference, AE, between the pho- 
tosensitive units, (ii) on the Forster coupling Vp, and 
(iii) on the Forster reorganization energy A^ . The large 
energy separation, AE, of the energies of the nearby pho- 
tosensitive molecules leads to a decrease of the Marcus 
rates and, thus, to the suppression of the energy transfer. 
In our case, the energy distance between the BPEA (A 
= 450 nm) and BDPY (A = 513 nm) molecules is about 
340 meV, whereas the energy separation of the BDPY 
chromophore and porphyrin (A = 620 nm) is of order of 
415 meV. This energy gap can be partially compensated 
by the large reorganization energy, A p , which reflects the 
significant fluctuations of the relative positions of the en- 
ergy levels. Here, the environment plays a positive role 
assisting the efficient and fast energy transfer between 
chromophores (see also Ref. [HI [JO])- The timescales for 
the energy and electron transfers should be shorter than 
the radiation leakage time and the quenching time, oth- 
erwise the energy of the photons absorbed by the system 
will be lost. 

At very low temperatures (e.g., liquid nitrogen, T = 
77 K), the fluctuations in the positions of the energy 
levels of the chromophores are frozen, so that the light- 
to-electricity conversion requires sufficiently large values 
of the reorganization energy, Ap > 180 meV (see Figs. 6a 
and 7a). Note that for Vp > 25 meV the linear-chain 
arrangement system has an optimal performance at Ap — 
225 meV for both frequency ranges. This extreme low 
temperature case is only shown for comparison with the 
higher temperatures cases. 

At room temperature (T = 298 K) and at strong 
enough Forster coupling, Vf > 25 meV, the blue and 
green spectral peaks demonstrate similar behaviors as 
functions of Ap (see Figs. 6b and 7b). Here, the quantum 
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yield begins to grow when the reorganization energy ex- 
ceeds ~ 80 meV, reaching finally 90% at Ap > 150 meV. 
These numbers are determined by the parameters of the 
antenna-triad complex, and, especially, by the radiation 
leakage time r ra d, of the excited porphyrin state P*, es- 
timated above as T ra( j ~ 0.4 ns. 

At high temperatures (see Figs. 6c and 7c, plotted for 
T = 500 K) the facilitating effect of the environment 
increases, and the efficient energy transfer starts at the 
lower reorganization energies, Ap > 75 meV. This value 
of Ap is comparable with the reorganization energy for 
the energy transfer in the B850 light-harvesting complex, 
where A F ~ 27 meV [Tg] . 

At the smaller value of the Forster coupling, Vp = 
1 meV, the conversion of the blue light (A = 450 nm) to 
electricity is significantly suppressed (see Fig. 6) , whereas 
for green light (A = 513 nm) the dependence of the quan- 
tum yield on Aj? are shifted to higher reorganization ener- 
gies (Fig. 7), compared to the case of the larger couplings, 
V F = 25, 50 meV. 

It should be noted that, to cover a broader range of the 
spectrum of light with a fixed number of antenna chro- 
mophores, the resonance energies of the light-harvesting 
complexes should be very well separated. However, in 
this case the energy transfer between the antenna chro- 
mophores would be quite slow, since this transfer is gov- 
erned by the rates corresponding to the inverted regions 
of the Marcus parabola. Then, the dissipation comes into 
play, and the energy of the absorbed photons is lost on 
its way from the antennas to the reaction center. The en- 
ergy transfer rates and, thus, the efficiency of the system 
can be maximized in the case when the energy distance, 
AEi = E i+ i — Ei, between the nearby light-harvesting 
complexes (labeled by indices i + 1 and i, with energies 
Ei+i and Ei) is equal to the corresponding reorganization 
energy, A F . That is, AEi = E i+1 - E t = A F . 



VI. CONCLUSIONS 

In this paper we have studied theoretical aspects of 
the operation of an artificial reaction center (a ferrocene- 
porphyrin-fullerene molecular triad) coupled to the com- 
plex of four light-harvesting molecules. We have ana- 
lyzed two configurations of the antenna complex: (a) 
a star-shaped configuration, where each light-harvesting 
molecule is able to transfer energy directly to the 
centrally-located reaction center, and (b) a case where 
the antenna molecules form a linear chain, which gradu- 
ally transfers excitations from the high-energy antenna 
located in the far end, to the antenna chromophore 
with the lowest energy. The last antenna chromophore 
in the chain is energetically connected to the reaction 
center (RC). To be specific, we have considered the 
case when the antenna complex is comprised of two 
molecules of bis (phcnylethynyl) anthracene (BPEA), ab- 
sorbing blue photons (A = 450 nm), and two molecules 
of borondipyrromethene (BDPY), having an absorption 



maximum in the green region (A = 513 nm). We have 
shown that the configuration with a linear arrangement 
of the antenna chromophores (configuration (b)) is able 
to convert blue and green photons to electricity with a 
quantum yield of order of ^30% (over a wide range of 
wavelengths), whereas the energy of the red photons, ab- 
sorbed by the molecular triad itself (A = 620 nm), is 
converted to a current with a quantum yield reaching the 
value of 90%. We have investigated dependencies of the 
quantum yield on the Forster reorganization energy as 
well as on the Forster coupling constants between chro- 
mophores and have shown that the environment plays 
a significant role in facilitating the antenna-RC energy 
transfer, thus, improving the light-harvesting function of 
the system. Overall, the configuration (b) is more effi- 
cient than (a) in transferring energy to the reaction cen- 
ter. 

We emphasize that the artificial photosystem an- 
alyzed in this work can be implemented with real 
light-harvesting components, such as porphyrin and 
BPEA/BDPY molecules. The excitonic (Forster) cou- 
pling strongly depends on the mutual distances and the 
orientations of the chromophores. Similar to the wheel- 
shaped antenna-reaction center complex implemented in 
Ref. [5|, the components of the photosystem can be 
placed at distances of the order 10 A, which allows for a 
sufficiently strong Forster coupling between the antenna 
chromophores and the reaction center. At the same time, 
the chromophores comprising the light-harvesting com- 
plex retain their individual molecular features. The re- 
organization energy, another controlling parameter for 
energy transfer, is varied for the system under study. 
Namely, we numerically calculate both the light-induced 
electron current and the quantum yield as functions of 
the reorganization energy. This allows us to determine 
the value of the reorganization energy at which the sys- 
tem works with maximum optimal efficiency. 
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Appendix A: Hamiltonian 

Here we describe the methods used in our work. We 
characterize the electrons in the states i (= D, P, P*, 
Anl, Anl*, An2, An2*, An3, An3*, An4, An4*, A) by the 
Fermi operators af and with the electron population 
operator rii = a^ai. Each electron state can be occupied 
by a single electron, as the spin degrees of freedom are ne- 
glected. Electrons in the leads (electrodes) are described 
by the Fermi operators dt , dk, a , where a = L, R; and k 
is an additional parameter which has the meaning of a 
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wave vector in condensed matter physics. The number 
of electrons in the leads is determined by the operator 
^2 k Nk a , with N ka = d^ a d ka . The total Hamiltonian 
of the system is complicated. It includes the terms de- 
scribed below. 



1. Eigenenergies and Coulomb interactions 

This part of the Hamiltonian involved the eigenener- 
gies of the electron states (i — D, P, P*, Anl, Anl*, 
An2, An2*, An3, An3*, An4, An4*, A) and the Coulomb 
interactions between the electron states. 

H a = EjUj + upnpnp, + u DP (l - n D ) (1 - n P - 

i 

n P .) - u DA {\ - n D )n A - u PA {\ - n P - n Pif )n A . 

(Al) 

The symbols up,up,p,up, Al up A represent the electro- 
static interactions between the electron sites. We have 
assumed that the empty donor state D (with np> = 0) as 
well as the empty photosensitive group (up + np* =0) 
have positive charges. Therefore, Up>p > because both 
D and P are positively charged and thus repulsive. The 
acceptor state A becomes negatively charged when it 
is occupied by an electron and thus —Uda < and 
— Up A < 0. This attraction occurs when the acceptor 
A is occupied (n A = 1) and the D and P states are both 
empty. Also, the acceptor state A is neutral when it is 
empty. 



2. Forster couplings 

We consider the energy transfer between the reaction 
center (P) and the antenna complexes, and also among 
the antenna complexes by introducing Forster coupling 
terms, 



-^Forster = - J~] Vkl a ( » ajj". a k + H.c. 



(A2) 



kl 



where, the pair {k,l} = {P, Anl}, {P, An2}, {P, An3}, 
{P, An4}, {Anl, An2}, {An2, An3}, {An3, An4}. Here, 
Vf determines the strength of the Forster coupling and 
is proportional to the product of the dipole moments, er\ 
and er k , and inversely proportional to the cubic power of 
the separating distance, R, between the chromophores, 



V F = 



27T£ e R 3 



(A3) 



3. Tunneling couplings to the leads 

The electron tunnelling from the left lead to the donor 
state and from the acceptor state to the right lead are 



both given by the Hamiltonian, 

Htr = - TkL a D C kL - TkR C kR a A + H ' C - » ( A4 ) 

k k 

where c ka , c ka are the electron creation and annihilation 
operators, and a is the index for the leads. The Hamil- 
tonian of the leads is given by 



Hlr — ^ £ a n a 



with 



— c^ a c ka . 



4. Thermal tunneling 

Activated by thermal fluctuations, electrons can tun- 
nel between the electron sites. Here, H tU n, given by the 
following expression, 



#tun = - ^ A M / a+ a v + H.c. , 



(A5) 



which accounts for the thermal tunneling effects. Here 
Ay is the strength of the tunnelling coupling and the 
{1,1'} indices refer to the pairs: {D, P}, {D, P*}, {A, P}, 
{A, P*}. 



5. Light-induced excitations 

This part of the Hamiltonian accounts for the inter- 
action of light with the molecular triad and the antenna 
complexes. Under the rotating-wave approximation, the 
light-induced excitation processes can be described as 



H 



Light 



Y,Fe iuJot a+a k . + H.c. (A6) 
k 



where k = P, Anl, An2, An3, An4 and the field amplitude 

F £cxt d kk * , 
where, d kk * is the dipole moment. 



6. Coupling to a radiation heat bath and an Ohmic 

bath 



Coupling the system to a radiation heat bath causes 
radiation leakage from the excited states. The following 
Hamiltonian accounts for this radiation leakage 



H Q = - ^2 Qoo'ata a , + H.c. 



(A7) 



where, {a, <r'} denotes the pairs of sites {D, P*}, {A, 
P}, {P, P*}, {Anl, A nl*}, {An2, An2*}, {An3, An3*}, 
{An4, An4*}. 

The operators for the radiation bath, 



(A8) 
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are proportional to the projection of the fluctuating elec- 
tromagnetic field, £rad) along the direction of the corre- 
sponding dipole moment, d a(J i — e x aa >. 

The excited state of the photosensitive part of the 
molecular triad can be quenched by the electrode. 
Namely, lose the excitation energy when interacting with 
the electrodes. We introduce i? quenc h to account for this 
energy- loss or quenching processes., 



H, 



quench 



-Qi at a v + H.c. , 



(A9) 



where, Qi is the variable of the Ohmic bath and {I, I'} — 
{P, P*}. 



7. Interaction with the environment 

We have taken into account the effects of a dissipative 
environment by the well-known system-reservoir model 

HHIT31 na na 



2m a 



-Y 

2 ^ 



(A10) 



where Xj , Pj are the position and momentum of the jth 
oscillator with effective masses rrij and frequencies ujj. 
Here, Xji is a measure of the strength of the coupling be- 
tween the electron subsystem and the environment. We 
characterize the phonon modes of the bath by the spec- 
tral functions J^^lo), defined by 



J - 5{u - uj) . (All) 

3 



The spectral function Ju> is related to the reorganiza- 
tion energy Aw for the i — > i' transition, by the following 
equation: 

A,;, = / " *Lj„( U ) = £ ™>^ -*>«)\ (A12) 

3 



8. Unitary transformation 

To remove the electron population operators of the 
electron subsystem from H env , we use the unitary trans- 
foration U — Yli Ui , where 



Ui = exp 



E 



(A13) 



The results of this unitary transformation are: 

(i) When U is operated on an arbitrary function <&(xj), 

a shift of the oscillator positions xj is produced. 



(A14) 



(ii) Another result of this unitary transformation is that 
all the transition operators acquire fluctuating factors. 

tfalavU = e*»'a\a v . (A15) 
The stochastic phase operators £u> are given by 



tw = & - &' with 6 = i E 



PjXji . 



(A16) 



Appendix B: Master equations 

The system under study can be characterized by the 
256 eigenstates of H a . We expressed all the operators 
described in previous section in the terms of the density 
operators p^ = \p)(v\. To derive the time evolution 
of the diagonal elements p MM = p^ of the density matrix 
(Pnu)i we write the Heisenberg equation for the operators, 
with the subsequent averaging (p^) over the environment 
fluctuations. Finally, we obtain the master equation for 
the density matrix of the system [TTJ [T21 [TSl El [55] , 



(Bl) 



where, 7^ is the total relaxation matrix, which is the 
sum of six types of relaxation rates: 

= 7^ + fe^ rster + fcJSf + ^ ht + Kl d + k T" lch ( B2 ) 

These relaxation rates will be described right below. 



1. Electron tunneling rates between the leads and 
the molecular triad 



The first term of Eq.(B2|, 7* r „, represents the relax- 
ation rates due the couplings of the triad to the L and R 
electron reservoirs: 

1% = t l {\a Dl ^\ 2 [l - + \aD;vn\*fL(u^)} 

+T R {\a A .^\ 2 [l - /r(w„ m )] + |o^| 2 /r(w pw )} 03) 

where the r a (a = L, R) are the resonant tunneling rates. 
Here, the electron reservoirs have been characterized by 
the Fermi distributions / Q (w), 



fa{Eka) 



exp 



E 



k B a 



Pa 



T 



-i -1 



(B4) 



with the temperature T (fee = !■,%= 1). The electro- 
chemical potentials pl and pn are the controlling fac- 
tors of the electron transition rates from the left lead to 
the donor state, and from the acceptor state to the right 
lead. 
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Forster relaxation rates 



5. Relaxation rates due to radiation leakage 



Here fc^° rstGI accounts for the excitation transfer rates 
from the antenna complexes to the reaction center, and 
also among the antenna complexes. The excitation tran- 
sition rates via the Forster mechanism is given by 



Forster 



1/ + + \ 1 21 T + A-kl) 2 

\(afai*al,a k ) uli \ \ exp 4A f — 



(B5) 

where Vki is the resonant Forster relaxation rate and Afc; 
stands for reorganization energy. We denote, Vu = Vf 
and Am = Ap for any combinations of k and I. The non- 
resonant exponential term of the above expression arises 
due to the different energy gaps of the reaction center and 
the accessory antenna complexes. Moreover, the non- 
resonant exponential terms depend on the reorganization 
energy. 



Thermal tunneling rates 



The matrix element fc*" n of the Eq. (B2| are respon- 
sible for the relaxation processes arising from thermal 
tunneling. These are given by 



\(a+a a ') VIM \ 2 ] exp 



(oj^ + A aa >) 



- 1 



,(B6) 



where A is the resonant tunnelling rate, uj^ is the en- 
ergy difference between the states /i and v (acting as a 
thermodynamic gradient), and the reorganization energy 
A are the main guiding factors of the thermal tunneling 
rates 



c tun 



4. Light-induced excitation rates 



Neglecting the effects of the environment on the radi- 
ation transitions, fc^ d is given by 



„rad 



2n 



J2 l<W| 2 [|(a+a ff 0^| 2 + l(^<W| 2 ] 



(Sf)V(^) 



(B8) 



where, n and d aa i stand for the refraction index and the 
dipole moment, respectively. 



6. Lead-induced quenching rates of the excited 
states 



The last term of Eq. (B2), fc^° nch > describes the en- 
ergy loss due to the quenching of the excited state of the 
photosensitive part of the triad. 



<" ench = a p {\(a+a P ,)^\ 2 + \(a+a P ,) v ^ 



fcothf^) 


- 1 


. \2T J 





(B9) 



Appendix C: Current and efficiency 

1. Electron current 

For weak couplings, the electron flowing between the 
leads and the molecular triad is given by: 



dt 

We derive the equation of the current in terms of the 
density matrix elements. 

h = Fr^2 \ a A-,^\ 2 [l ~ fn{u VIJ )]{p v ) 

(XV 

~ \ a A;^\ 2 fR(^^^)(P^) ■ (CI) 



The contribution fcj]| ht to the total relaxation matrix 
due to light-induced excitation processes is 

2 



,.light 



= £1*1 

k 

x exp 



{|(a+a 



Afjcr* T 

{uJnv + cj + A^.) 2 



(a+a^^exp 



(u^ - cjq + A go .») 2 
4A„„.T 



(B7) 



This rate includes contributions from the following tran- 
sitions, P-^P*, Anl -» Anl*, An2 ->■ An2*, An3 ->■ An3*, 
and An4 — !• An4* . 



2. Absorbed energy 

The total amount of energy absorbed per unit time 
^photon by the molecular triad and antenna chromophores 
is 



^photon = £w |F| 2 J ^ ^ \( a * a <r)^\ 2 (P»- 



Pv 



x I exp 



exp 



(ujpv - A^* + UJ ) 2 

4A ffCT ,T 
{^nv — A aa * — wq) 2 



4A„„.T 



(C2) 



where a = P, Anl, An2, An3, An4. 
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Power-conversion efficiency The quantum yield is denned as 



The power-conversion efficiency of the system is the the , , ^ , , t . () 

ratio of the output -E ou tput and the input -Ei npu t energies, * = ~jy = — 77~ ■ 

(C4) 

' 1 'photon ' 



^ _ -E-output _ £pump _ Ir{HR — Hl) (C3) 
-^input ^-photon ^photon 
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